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Abstract 

Three  examples  of  reeent  investigations  eombining  CFD  analysis  and  wind  tunnel  tests  at  the  von  Karman  Institute 
for  Fluid  Dynamies  (VKI)  are  presented  to  illustrate  the  integration  proeess  between  experimental  tests  and 
numerieal  results.  For  eaeh  projeet,  the  requirements  and  the  methodology  were  different  but  the  final  goal  was  to 
obtain  an  optimized  aerodynamie  design.  The  first  projeet  presented  is  dealing  with  a  fluid- strueture  interaetion 
problem.  The  seeond  projeet  eoneems  the  aerodynamie  design  of  a  new  Belgian  polar  base  station  for  Antaretiea 
and  the  last  projeet  eoneems  the  design  of  an  ultra- streamlined  land  vehiele  for  the  next  Solar  World  Challenge. 

Key  words:  fluid  dynamies,  vehiele  aerodynamies,  wind  tunnel,  CFD,  validation,  polar  base,  fluid- stmeture 
interaetion 

Introduction 

Integration  of  CFD  and  experimental  results  in  aerodynamies  is  required  when  eaeh  method  alone  is  not  reliable, 
fast  or  aeeurate  enough  to  reaeh  the  final  objeetive  of  the  projeet.  The  advantage  of  eaeh  teehnique  is  then  used  in 
an  integrated  projeet.  CFD  analysis  ean  be  performed  for  example  in  a  short  period  of  time,  but  requires  a  validation 
when  the  design  falls  outside  the  range  of  standard  aerodynamie  design. 

Based  on  the  experienee  aequired  at  the  von  Karman  Institute  (VKI),  the  best  integration  is  obtained  when  the 
seientists  performing  the  wind  tunnel  tests  and  the  CFD  people  are  integrated  in  a  global  team  knowing  the  final 
objeetive  of  the  projeet  and  able  to  have  a  eritieal  view  on  the  results  gained  using  wind  tunnel  or  CFD  eode.  To 
illustrate  the  integration  proeess,  three  examples  of  reeent  investigations  eombining  CFD  and  wind  tunnel  tests  are 
presented.  For  eaeh  projeet,  the  requirements  and  the  methodology  are  different  but  the  final  goal  was  to  obtain  an 
optimized  aerodynamie  design. 

The  first  projeet  deals  with  the  improvement  of  the  performanee  under  wind  load  of  new  sun  proteetion  system  to  be 
plaeed  alongside  building.  Vibrations  on  existing  stmeture  were  observed  during  heavy  storms.  The  investigation 
has  started  using  wind  tunnel  tests.  The  wind  tunnel  tests  were  used  to  eharaeterize  the  fluid- stmeture  interaetion 
leading  to  the  vibrations.  A  new  geometrieal  design  of  the  blades  of  the  proteetion  system  has  been  performed  using 
CFD  before  a  final  validation  in  the  wind  tunnel.  The  seeond  projeet  eoneems  the  aerodynamie  design  of  the  new 
Belgian  polar  station  for  the  Antaretiea.  From  an  initial  shape,  the  building  has  been  modified  using  a  eombination 
of  CFD  and  wind  tunnel  tests  at  redueed  seale.  Snow  drift  and  snow  erosion  were  analyzed  experimentally.  The 
results  of  the  numerieal  simulations  were  used  to  understand  the  flow  pattern  around  the  building  and  to  obtain  the 
global  loads  on  the  building.  New  designs  were  tested  both  numerieally  and  experimentally.  Loeal  pressure 
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distributions  obtained  by  CFD  are  now  used  in  the  meehanieal  design  of  the  different  external  panels.  The  third  and 
final  projeet  deals  with  the  aerodynamie  design  of  an  ultra- streamlined  land  vehiele,  namely  a  new  solar  ear  for  the 
next  World  Solar  Challenge  that  will  take  plaee  in  Australia  in  Oetober  2007.  A  ehange  in  the  regulation  for  the 
new  raee  has  foreed  the  teams  to  ehange  the  overall  design  of  the  ears.  CFD  has  been  used  for  the  initial  design  and 
the  definition  of  the  ear  eoneept  but  the  final  geometrieal  eonfiguration  was  tested  and  optimized  with  the  help  of  the 
wind  tunnel  tests. 


Description  of  the  test  faciiities  and  CFD  code 


The  experimental  tests  for  the  three  investigations  are  performed  in  the  VKI  Low  Speed  wind  tunnel  LI.  The  LI 
wind  tunnel  is  a  elose-loop  wind  tunnel  built  during  the  fifties  and  eontinuously  upgraded.  The  two  eontrarotating 
fans  of  this  faeility  are  driven  by  a  variable  speed  DC  motor  of  590  kW,  allowing  a  eontinuous  variation  of  wind 
tunnel  veloeity.  The  eontraetion  ratio  of  the  settling  ehamber  is  4  with  a  typieal  turbulenee  level  of  0.3%.  Two  tests 
seetions  ean  be  adapted  for  this  wind  tunnel: 

•  The  open-jet  test  seetion  LI -A  (Figure  1)  has  a  free  jet  test  seetion  of  3  m  diameter  and  4.5  m  length  and 
was  designed  as  an  aeronautieal  wind  tunnel  and  is  used  for  testing  aireraft  models  with  spans  up  to  about  2 
m  and  2D  wings  at  Reynolds  number  up  to  4.10Vm.  This  open  test  seetion  is  also  used  for  tests  on  ground 
struetures  and,  in  the  eurrent  paper,  has  been  used  for  the  first  investigation  dealing  with  fluid- strueture 
interaetion. 

•  The  wind  engineering  faeility  Ll-B  (Figure  2)  eonsists  of  a  speeially  designed  test  seetion  of  2  m  high,  3  m 
wide  and  20  m  long  with  a  roughened  floor  to  allow  the  growth  of  a  turbulent  boundary  layer  similar  in 
nature  to  the  lower  part  of  the  atmospherie  neutral  boundary  layer.  It  is  therefore  an  ideal  tool  for 
mieroelimate  studies  around  buildings,  like  for  the  new  Belgian  polar  base  presented  as  seeond  example  in 
this  paper.  A  speeial  test  ehamber  is  also  provided  at  the  beginning  of  the  seetion  for  investigation  of 
surfaee  vehiele  or  boat  aerodynamies,  and  is  equipped  with  a  six-eomponent  floor  balanee  and  faeilities  for 
boundary  layer  suetion.  The  measurements  performed  on  the  solar  ear  were  made  using  this  balanee  at  the 
entranee  of  the  wind  tunnel,  meaning  that  the  flow  veloeity  was  uniform  when  reaehing  the  model. 


COWTRAROUTING  PROPEUERS 


MOTOR  CONTRAROTATING  PROPELLERS  DIFFUSER 


Figure  1  -  wind  tunnel  Ll-A 


Figure  2  -  wind  tunnel  Ll-B 


The  CFD  investigation  is  performed  using  the  eommereial  CFD  eode  FLUENT  from  the  eompany  ANSYS.  This 
well-know  eommereial  CFD  eode  allows  simulating  fluid  flow,  heat  transfer  and  ehemieal  reaetion.  The  eode 
solves  the  eonservation  equations  for  mass,  momentum,  energy  and  ehemieal  reaetion  using  a  finite  volume  method. 

Flow-structure  interaction 

The  first  projeet  deals  with  the  improvement  of  the  performanee  under  wind  load  of  new  sun  proteetion  system  to  be 
plaeed  alongside  building,  like  shown  in  Figure  3.  Vibrations  on  existing  strueture  were  observed  during  heavy 
storms,  and  our  elient  has  asked  us  to  modify  the  aerodynamie  shape  of  the  blades  to  eliminate  or  to  diminish  the 
fluid-strueture  interaetion. 
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The  investigation  started  with  wind  tunnel  tests  to  eharaeterize  the  fluid- strueture  interaetion,  and  more  preeisely  to 
get  the  eritieal  wind  veloeities  leading  to  vibration,  and  even  to  damages  and  destruetion  of  existing  wind  proteetion 
systems.  High  speed  image  reeording  during  the  tests  has  provided  the  deformation  mode  of  the  different  blades  as 
funetion  of  the  wind  veloeity.  Figure  4  shows  a  snapshot  taken  during  a  wind  tunnel  test.  The  deformation  of  a  high 
number  of  blades  ean  be  seen  on  this  figure.  By  analyzing  in  more  details  the  video  reeorded  during  the  wind  tunnel 
tests,  we  eonelude  that  the  vibrations  of  the  blades  are  more  related  to  torsion  than  bending.  Therefore,  from  the 
wind  tunnel  tests  we  eame  to  the  eonelusion  that  the  seeond  meehanieal  mode  is  the  one  exeited  by  the  flow  (wind). 


The  theoretieal  investigation  has  eombined  a  meehanieal  analysis  using  a  CAE  paekage  and  a  CFD  analysis.  The 
meehanieal  analysis  is  performed  to  eompute  the  natural  frequeneies  of  the  blades  for  different  geometrieal 
parameters  and  for  the  eomputation  of  the  deformation  of  a  blade  under  wind  load.  The  CFD  analysis  of  the  flow 
field  around  existing  blade  is  performed  to  get  the  vortex  shedding  frequeney  and  the  Strouhal  number  assoeiated  to 
the  unsteadiness  of  the  flow  field  [1].  Vibrations  of  the  blades  will  oeeur  when  the  frequeney  at  whieh  the  vortiees 
are  shed  from  the  blade  will  eorrespond  to  one  of  the  natural  frequeney  of  the  blade. 


Figure  5  shows  the  frequeneies  of  the  different  meehanieal  modes  as  funetion  of  the  blade  length  (distanee  between 
the  two  supports).  The  first  and  the  third  mode  eorrespond  to  a  bending  mode,  while  the  seeond  mode  eorresponds 
to  torsion.  As  predieted  by  the  theory  for  bending  and  torsion  modes,  a  linear  evolution  of  the  frequeneies  is 
obtained  as  funetion  of  the  square  of  the  inverse  of  the  distanee  between  the  supports  of  the  blades  (Figure  6).  A 
first  solution  to  eliminate  the  vibrations  observed  during  the  wind  tunnel  tests  would  be  to  shorten  the  distanee 
between  the  supports,  but  this  is  eommereially  not  aeeeptable  as  more  supports  are  needed  per  unit  length.  The 
solution  is  therefore  in  a  better  aerodynamie  shape  of  the  blades. 
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Figure  5  -  mechanical  analysis  -  blade  deformation 


Figure  6  -  Effect  of  the  length  on  the  vibration 
frequency 
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For  the  prediction  of  the  vortex  shedding  frequencies,  the  commercial  code  FLUENT  has  been  used  to  predict  the 
evolution  of  pressure  and  velocity  as  function  of  time.  Numerical  simulations  have  been  performed  using  a  second- 
order  implicit  time  stepping  for  the  transient  computations.  The  turbulence  is  modeled  using  a  R.A.N.S.  approach, 
with  a  k-8  realizable  turbulence  model  with  enhanced  wall  treatment.  Special  care  has  been  taken  in  the  grid 
generation  to  decrease  the  value  of  Y+  near  the  blade.  The  numerical  simulations  have  been  performed  on  a  2D 
cross-section  of  the  blade  to  stay  within  the  time  constraints  imposed  by  our  client. 

The  simulations  are  performed  for  the  same  wind  direction  with  respect  to  the  blades  as  for  the  wind  tunnel  tests. 
Contours  of  velocity  magnitude  at  two  different  time  step  are  shown  in  Figure  7  and  Figure  8,  showing  the  velocity 
fluctuations  as  function  of  time.  The  variations  of  drag  and  lift  are  periodic  and  the  vortex  shedding  frequency  is 
computed  for  different  wind  velocities.  The  results  are  shown  in  Figure  9.  A  Strouhal  number  has  been  associated  to 
the  periodic  fluctuations  and  a  constant  value  has  been  found  for  the  different  numerical  wind  velocities. 


For  a  wind  velocity  of  10  m/s,  the  vortex  shedding  frequency  found  numerically  is  equal  to  33  Hz,  which  is  close  to 
the  second  mode  (torsion)  of  the  blade  for  a  distance  between  the  supports  equal  to  1520  mm.  Vibrations  are 
therefore  predicted  to  start  at  a  wind  velocity  of  10  m/s.  This  result  has  been  validated  by  the  wind  tunnel  tests,  for 
which  vibrations  related  to  torsion,  were  observed  to  begin  at  a  wind  tunnel  velocity  of  10  m/s. 
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Figure  7  -  contours  of  velocity  Figure  8  -  contours  of  velocity  Figure  9  -  results  for  one  type  of  blade 
magnitude  at  t=1.125  s.  magnitude  at  t=1.13  s.  -  Strouhal  number 


The  validated  theoretical  approach  has  then  been  used  in  a  next  step  to  define  a  new  aerodynamic  design  of  the 
blades  of  the  sun  protection  system.  Ten  new  different  geometries  have  been  analyzed  theoretically,  both  on  a 
mechanical  point  of  view  using  the  CAE  code  and  on  the  fluid  dynamics  aspects  using  the  CFD  code. 

The  natural  vibration  frequencies  of  the  ten  new  blades  were  very  close  to  each  other  for  the  three  first  modes,  and 
therefore  the  best  shape  of  blade  had  to  be  found  using  unsteady  numerical  simulation  of  the  flow  field  around  the 
blade.  In  order  to  reduce  the  amplitude  of  the  fluctuations  in  lift  obtained  by  CFD,  the  aerodynamic  shape  of  the 
blade  was  modified.  In  the  final  design,  the  velocity  increase  on  the  upper  part  of  the  blade  has  been  considerably 
decreased,  leading  to  lower  pressure  fluctuations  on  the  surfaces  of  the  blade.  It  was  believed  that  this  would  result 
in  lower  amplitude  for  vortex  shedding  and  that  the  mechanical  mode  of  the  blade  would  not  be  excited  by  the  flow 
field. 

The  new  blade  with  the  newly  defined  aerodynamic  shape  has  been  constructed  and  a  final  validation  of  the  new 
design  was  performed  in  our  El -A  wind  tunnel.  Very  good  results,  meaning  no  vibrations  at  low  speed,  have  been 
obtained  with  the  proposed  blade  geometry. 

The  methodology  for  this  investigation  showing  the  integration  of  experimental  wind  tunnel  tests  and  numerical 
results  is  summarized  in  Figure  10.  The  investigation  started  with  a  wind  tunnel  test  reproducing  a  real  test 
condition  and  ended  with  the  same  wind  tunnel  test  validating  the  new  geometry.  The  investigation  would  not  have 
been  successful  without  the  combination  of  wind  tunnel  tests  and  CFD  simulation. 
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Figure  10  -  integration  of  wind  tunnel  tests  and  CFD  simulation  -  fluid/structure  interaction  study 


Aerodynamic  design  of  a  new  poiar  station 

The  second  project  concerns  the  aerodynamic  design  of  the  new  Belgian  polar  station  for  the  Antarctica.  This  new 
Belgian  Research  Station  will  be  constructed  in  the  Sor  Rondane  Mountains,  Antartica,  during  the  austral  summer 
2007-2008.  The  station  will  have  a  hybrid  design,  with  the  main  building  above  ground-level  (Figure  11)  and 
anchored  onto  a  snow-free  rock  area  (existing  ridge).  The  presence  of  the  most  intense  winds  on  Earth,  so  called 
katabatic  wind  leading  to  wind  speed  in  excess  of  65  m/s,  transport  large  quantities  of  snow  that  can  bury  the 
building  very  quickly.  Two  important  aspects  are  therefore  to  be  considered  in  the  aerodynamic  design  of  an 
Antarctic  building:  snowdrift  and  wind  loading. 


Figure  11  -  artistic  impression  of  the  future  base  Figure  12  -  snow  erosion  in  wind  tunnel  behind  the 

building 


Today,  almost  all  modem  Antartic  stations  have  undergone  aerodynamic  studies  at  different  stages  of  design. 
References  on  recent  aerodynamic  studies  of  new  Antarctic  stations  (Neumayer  III,  Halley  VI)  can  be  found  in  the 
article  of  J.  Sanz  [2].  For  the  new  Belgian  polar  station,  an  aerodynamic  study  has  been  included  from  the  beginning 
of  the  building  design  process. 

The  aerodynamic  design  of  the  new  Belgian  station  at  VKI  is  divided  in  two  phases: 

•  A  conceptual  design  phase 

•  A  detailed  design  phase 
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In  the  conceptual  phase,  we  have  analyzed  the  ability  of  different  buildings  to  cope  with  intense  drifting  snow.  As 
mentioned  previously,  many  Antarctica  stations  have  had  serious  problems  with  snow  drifting  with  eventual  buried 
structures  that  require  costly  maintenance  to  prevent  inaccessibility  and  structural  damage.  The  use  of  CFD  for  the 
evaluation  of  snow  drifting  and  snow  deposition  was  rejected  from  the  beginning  because  it  would  have  involved  the 
development  of  new  models  that  would  require  validation  using  wind  tunnel  tests.  The  timeframe  of  the  conceptual 
phase  did  not  allow  performing  this  development.  Therefore,  the  evaluation  of  the  snow  drifting  and  snow 
deposition  around  different  building-block  shapes  has  been  made  exclusively  using  wind  tunnel  tests.  However, 
snowdrift  modeling  in  wind  tunnel  is  still  a  difficult  task  due  to  the  inherent  impossibility  to  maintain  similarity  of 
all  the  driving  forces  and  the  experimental  study  of  snow  drifting  is  also  very  time  consuming.  As  a  result,  it  was 
not  feasible  to  use  this  technique  as  many  building  configurations  had  to  be  tested  in  a  very  limited  amount  of  time. 
Instead,  the  sand  erosion  technique  proved  to  be  a  very  efficient  tool  to  evaluate  the  wind  conditions  at  ground  level 
and  was  used  to  point  out  areas  of  severe  wind  comfort  and  at  the  same  time  identify  snow  accumulation  and  erosion 
regions.  Of  course,  the  sand  erosion  technique  is  a  surface  visualization  and  does  not  provide  information  about  the 
development  rate  of  the  buildup  height  (snow  accumulation).  However,  in  the  methodology  developed  and  based  on 
the  comparison  of  erosion  speed  contour  maps  with  and  without  the  building,  leading  to  an  amplification  factor  A, 
taken  as  a  measure  of  the  action  of  the  building  on  the  flow,  has  allow  VKI  to  find  regions  of  snow  accumulation 
and  snow  erosion.  More  information  on  the  technique  can  be  found  in  the  article  of  J.  Sanz  [2]. 

Figure  12  shows  a  typical  result  of  an  erosion  tests  with  the  polar  station  installed  on  the  ridge  topography  at  scale 
1:100  in  the  Ll-B  atmospheric  boundary  layer  wind  tunnel.  The  sand  erosion  produced  by  the  vortices  behind  the 
building,  where  the  garages  of  the  base  will  be  placed,  can  clearly  been  seen  on  this  figure. 

The  conceptual  design  phase  has  lead  to  the  evaluation  of  six  building  block  concepts  that  were  the  result  of  a 
preliminary  trade-off  between  different  design  parameters,  with  aerodynamics,  energy  demands  and  efficiency 
playing  the  leading  role.  At  the  end  of  this  phase,  the  most  suitable  basic  design  was  obtained  that  has  been 
finalized  in  the  detailed  design  phase. 


Figure  13  -  numerical  simulation  showing  the 
vortices  behind  the  building 


Figure  14  -  ‘‘Flow  visualization”  using  CFD 


In  the  detailed  design  phase,  the  height  and  the  position  of  the  building  on  the  ridge  has  been  defined.  Aerodynamic 
wind  loading  (drag,  lift  and  lateral  forces)  on  the  building  has  been  evaluated  using  CFD  for  different  building-ridge 
configuration.  A  typical  CFD  result  is  shown  in  Figure  14,  showing  the  pressure  distribution  on  the  building 
together  with  the  velocity  field  and  streamlines.  The  results  of  the  numerical  simulations  have  been  used  to 
understand  the  flow  pattern  around  the  building  (Figure  13)  and  to  obtain  a  first  estimate  of  the  global  loads  on  the 
building.  The  CFD  simulations  have  been  performed  for  one  wind  direction,  namely  perpendicular  to  the  building 
and  the  ridge. 

In  the  meantime,  a  scaled  model  of  the  final  design  instrumented  with  pressure  tabs  has  been  constructed  for  a  final 
validation.  The  final  building  is  shown  in  Figure  15  with  159  1mm  diameter  pressure  taps  connected  for  four 
scanning  valves.  The  pressure  measured  by  the  four  pressure  transducers  are  scaled  by  the  freestream  velocity 
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measured  at  a  referenee  position  and  maps  of  pressure  eoeffieient  have  been  obtained  on  the  different  surfaees  of  the 
building  for  different  wind  direetions. 

Loeal  pressure  distribution  obtained  by  CFD,  like  the  ones  shown  in  Figure  16  eould  then  be  validated  by  the  wind 
tunnel  experiments  and  the  loeal  data  from  the  CFD  eode  ean  then  be  used  for  the  meehanieal  design  of  the  external 
panels  of  the  polar  base. 


Figure  15  -  final  model  -  instrumented  with  pressure 
tabs 


Figure  16  -  CFD  results  showing  the  local  pressure 
distribution 


The  methodology  used  in  this  investigation  showing  the  integration  of  experimental  and  numerieal  results  leading  to 
useful  results  for  the  final  design  of  the  polar  base  is  summarized  in  Figure  17.  Wind  tunnel  tests  have  been  used  at 
the  beginning  of  the  study  for  predesign  studies  and  at  the  end  for  the  validation  of  the  final  design.  In  the  mean 
time,  CFD  has  been  used  to  define  the  position  of  the  base  on  the  ridge  and  to  understand  the  flow  field  behind  the 
building. 


WIND  TUNNEL  TESTS 


CFD  SIMULATION 


Predesign  of  the  building  - 
conceptual  design  phase 

Flow  visualization  using  CFD 

Database  on  snow  accumulation 
and  snow  erosion  for  typical 
building  shapes 


Determination  of  force  and 
moments  on  the  building 


Parametric  investigation  of 
the  position  of  the  station 


Wind  tunnel  tests  on  final 
geometry 


CFD  simulation  on  final  geometry 


Determination  of  local  loads 


Figure  17  -  methodology  -  design  of  a  new  polar  base 
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Aerodynamic  design  of  an  uitra-streamiined  iand  vehicie 


The  World  Solar  Challenge  is  organized  every  two  years  and  is  a  3000  km  long  raee  for  ears  powered  by  the 
sunlight.  In  2005,  the  winner  eompleted  the  journey  with  an  average  speed  in  exeess  of  100  km/h.  Students  from  the 
Engineering  Sehool  of  Leuven  (Belgium)  have  been  the  first  Belgian  team  to  partieipate  in  2005  with  a  solar  ear 
named  Umiear.  They  finished  the  raee  at  the  eleventh  position.  The  initial  aerodynamie  design  and  the  validation 
using  wind  tunnel  tests  were  performed  by  the  students  with  the  help  of  VKI  [3]. 


Faeed  by  the  possibility  that  speed  limits  may  be  imposed  for  the  2007  ehallenge;  a  new  generation  of  solar  ear  was 
ereated.  A  solar  ear  whieh,  with  little  modifieation,  eould  be  the  basis  for  a  praetieal  proposition  for  sustainable 
transport.  The  new  rules  inelude  a  limit  of  6  sq  meters  for  the  solar  eolleetors,  an  easy  driver  aeeess  and  an  upright 
seating  position.  These  rules  are  affeeting  the  external  aerodynamie  shape  of  the  solar  ear  and,  to  partieipate  to  the 
2007  World  Solar  Challenge,  students  from  the  Engineering  Sehool  Groep  T,  together  with  the  VKI,  have  studied 
and  defined  the  aerodynamie  shape  of  a  new  vehiele. 


Figure  18  -  computational  domain  +  surface  grid 


Figure  19  -  contours  of  pressure  -  CFD  results 


The  design  of  the  new  ear  was  started  using  the  well-known  referenee  book  for  the  aerodynamie  design  of  ultra- 
streamlined  land  vehiele  [4].  Numerieal  simulations  using  the  eommereial  CFD  eode  have  been  performed  from  the 
beginning  of  the  design,  first  in  2D  and  then  in  3D,  for  the  studies  of  different  aerodynamie  shapes  and  the  definition 
of  an  initial  eoneept.  In  this  first  phase  of  aerodynamie  design,  the  exaet  value  of  the  drag  was  not  the  major 
eoneem,  but  we  were  more  interested  in  trends  to  take  quiekly  deeision  on  for  example  the  external  shape  of  the  ear, 
on  the  number  and  position  of  the  wheels,  as  a  time  period  of  three  months  was  devoted  for  the  definition  of  the 
ear’s  eoneept.  The  eoneept  of  the  ear  at  the  end  of  this  first  phase  is  shown  in  Figure  18.  The  ear  will  have  two 
front  wheels  and  will  be  powered  by  the  baek  wheel  below  the  driver. 

After  this  first  phase,  a  seeond  phase  was  started  with  more  detailed  analysis  of  the  flow  field  around  the  ear  and  the 
influenee  of  loeal  geometrieal  modifieation  on  the  aerodynamie  of  the  ear,  like  for  example  the  thiekness  of  the 
wheel  system,  the  shape  of  the  eanopy,  the  position  of  the  eanopy  on  the  ear,  the  overall  shape,  thiekness  and 
eurvature  of  the  main  plate  supporting  the  solar  eells. 

A  typieal  result  of  the  seeond  phase  is  shown  in  Figure  18  (surfaee  grid)  and  Figure  19  (eontours  of  pressure  on  the 
ear).  A  total  maximum  number  of  5.550.000  eells  are  used  to  diseretise  the  eomputational  domain.  A  symmetry 
plane  is  used  to  limit  the  number  of  eells. 

During  this  seeond  phase,  different  numerieal  simulations  have  been  performed  by  ehanging  the  grid  distribution 
elose  to  the  body  of  the  ear  to  modify  the  eell  size  in  the  boundary  layer  and  also  by  using  different  turbulenee 
models.  Simulations  have  been  performed  with  the  k-e  realizable  turbulenee  model  and  with  the  Spalart-Allmaras 
turbulenee  model.  A  seatter  in  the  numerieal  aerodynamie  data  has  been  obtained.  For  example,  the  values  for  the 
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drag  area  coefficient  (Cd.A  in  m^)  obtained  by  numerical  simulation  were  between  0.12  and  0.15,  depending  on  the 
grid  distribution  and  the  turbulence  model  used.  This  scatter  is  quite  high  and  this  value  of  drag  coefficient  cannot 
be  used  to  predict  the  aerodynamic  performance  of  the  car,  and  by  this  way  compute  the  performance  of  the  solar 
car.  Tests  in  the  wind  tunnel  were  needed  to  get  more  accurate  aerodynamic  coefficients. 

At  the  end  of  the  CFDD  analysis,  a  first  overall  shape  was  defined  but  the  following  questions  were  remaining  and 
had  to  be  answered  using  wind  tunnel  tests: 

•  What  is  the  exact  value  of  the  Cd.A?  An  accurate  value  is  mandatory  to  predict  the  performance  of  the 
solar  car  and  to  develop  the  strategy  for  the  race. 

•  How  to  improve  slightly  the  geometry  of  the  canopy  to  reduce  the  overall  drag? 

•  What  is  the  best  geometry  for  the  wheel  protection  system  (thickness  -  length  of  the  fairing)  ? 

•  What  is  the  influence  of  side-wind  on  the  drag  and  side  force? 

•  What  would  be  the  influence  of  dust  or  dirt  on  the  surface  of  the  solar  cells  on  the  aerodynamic  drag? 

•  How  are  the  vertical  force,  pitching  moment  and  drag  changing  by  modifying  the  angle  of  attack  of  the  car 
by  a  few  degrees?  Modifications  during  the  race  will  be  possible  by  having  some  flexibility  on  the  shock 
absorber,  but  the  stroke  of  the  shock  absorbers  has  to  be  defined  in  this  second  phase. 

Wind  tunnel  tests  had  to  be  performed  because  it  was  the  only  way  to  get  in  a  short  timeframe  answers  to  these 
important  questions.  Once  the  different  models  for  the  wind  tunnel  were  built,  all  the  tests  have  been  performed  in 
one  week. 

The  wind  tunnel  tests  have  been  performed  using  two  different  scaled  models  at  a  scale  1/3,  meaning  that  by 
performing  tests  at  the  full  speed  of  the  wind  tunnel  (60  m/s),  the  Reynolds  number  could  be  reproduced  in  the  wind 
tunnel  for  a  car  driving  at  20  m/s. 

In  ultra- streamlined  vehicle,  most  of  the  drag  is  produced  by  skin  friction,  and  therefore  special  care  was  taken  to 
have  external  surfaces  as  smooth  as  possible.  A  front  view  of  the  solar  car  in  the  wind  tunnel  Ll-B  mounted  on  the 
six-component  floor  balance  is  shown  in  Figure  20.  Figure  21  shows  a  typical  result  obtained  with  oil  visualization 
performed  on  the  canopy  for  the  determination  of  the  detachment  position  of  the  flow  field.  This  result  has  been 
compared  with  the  prediction  of  the  flow  detachment  on  the  canopy  using  CFD,  and  a  very  good  agreement  was 
found  between  the  experimental  and  numerical  results. 

It  is  worth  to  mention  that  the  drag  area  coefficient  Cd.A  measured  during  the  experimental  tests  was  lower  than  the 
one  predicted  by  the  numerical  simulation  (reduction  by  more  than  10%).  This  is  one  of  the  very  good  news 
obtained  during  the  wind  tunnel  tests.  It  has  been  attributed  to  a  low  performance  of  the  CFD  code,  because  the 
numerical  model  was  not  tuned  for  the  prediction  of  the  transition  from  laminar  to  turbulent  of  the  boundary  layer. 
The  repeatability  of  the  wind  tunnel  tests  was  good,  as  the  scatter  in  the  Cd.A  coefficient  for  five  different  tests  was 
below  1%. 


Figure  20  -  ultra-streamlined  land  vehicle  in  wind 
tunnel  on  a  six-component  floor  balance 


Figure  21  -  Oil  visualization  showing  the 
detachment  position  on  the  canopy 
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The  methodology  of  the  design  proeess  with  the  integration  of  experimental  and  numerieal  results  leading  to  a  new 
ultra-streamlined  land  vehiele  is  summarized  in  Figure  22.  The  final  validation  of  the  aerodynamie  design  will  be 
made  during  the  raee,  for  whieh  the  expeetations  from  the  eurrent  team  are  quite  high.  A  position  on  the  podium  is 
expeeted. 


Figure  22  -  methodology  combining  CFD  and  wind  tunnel  tests  -  design  of  ultra-streamlined  vehicle 


Conclusions 

Using  three  different  examples  of  low-speed  aerodynamie  design,  eombining  wind  tunnel  experiments  and  CFD 
eomputation,  we  have  shown  how  integration  of  both  methods  ean  lead  to  a  better  aerodynamie  design.  For  eaeh 
projeet,  the  requirements  and  the  methodology  were  different,  but  the  final  goal  was  to  obtain  an  optimized 
aerodynamie  design. 

We  believe  that  the  best  integration  is  obtained  when  the  seientists  performing  the  wind  tunnel  tests  and  the  CFD 
people  are  integrated  in  a  global  team  knowing  the  final  objeetive  of  the  projeet  and  able  to  have  a  eritieal  view  on 
the  results  gained  using  eaeh  method.  Communieation  between  wind  tunnel  speeialists  and  CFD  speeialists  is 
therefore  erueial  throughout  the  study. 

The  three  examples  have  also  shown  that  wind  tunnel  tests  are  still  a  very  useful  tool  neeessary  for  aerodynamie 
design. 
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Introduction 


■  of  CFD  and  Experiments  is  when  is 

to  reach  the  final  objective  (with  correct  reliability,  delay  or  accuracy). 

■  The  of  each  method  is  then  in  an  integrated  project. 

■  The  is  obtained  when  Experimental  and  Numerical 

each  others. 


■  Illustration  through  three  examples: 

> 

> 

> 

■  For  each  project,  the  requirements  and  the  methodology  are  different  but  the  final 
goal  was  to  obtain  an  optimized  aerodynamic  design. 
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Description  of  the  test  facilities 

and  CFD  code 


■  Experimental  tests  are  performed  in  the  VKI  Low  Speed  Wind  Tunnel  L1 . 


CONTRAROTATING  PROPELLERS 


DIFFUSER 


SETTLING  chamber! 


HONEYCDMEl 


2.5m 


OREQANGUUR  TEST  SECTION  ^ 
L-IB  (WIND  ENGINEERING)  ^ 


:  (WIND  ENGINEERING) 

■  ■ 


CONTRACT  ON 


Collector 


FREE  JET  CIRCJLAR 
TEST  SECTION  L-lA 


L1-A  wind  tunnel  with  a  free  jet  test  section  L1-B  wind  tunnel  -  wind  engineering  facility 


>  Flow-structure  interaction  >  Test  of  the  Belgian  polar  station 

>  Test  on  the  ultra-streamlined  land  vehicle 


■  The  CFD  investigation  is  performed  using  FLUENT  from  the  ANSYS  with 
RANS  or  URANS  approaches. 
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system  to  be  placed  s  i  1 1  ud  :  bud  li  i  i. 
were  observed  heavy  st  r  ns. 

Modify  the  blade  shape  to  Is/diminish  the  1  j  I  s 


Flow  structure  interaction 
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Flow  structure  interaction 


■  The  investigation  has  tests. 

■  recording  has  provided  the  of  the  different 

blades. 
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Flow  structure  interaction 

■  Mechanical  investigations  using  a  Computer-Aided  Engineering  (CAE)  package 
(  )  has  the  natural  of  the  blades,  and  the 


Length  | 

#  supports 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

Mode  #5 

520mm 

2 

137.0Hz 

188.3Hz 

226.7Hz 

274.0Hz 

288.4Hz 

1000mm 

2 

50.2Hz 

62.7Hz 

153.4Hz 

155.3Hz 

186.5Hz 

1520mm 

2 

22.7Hz 

33.1Hz 

78.6Hz 

91.0Hz 

98.4Hz 

1800mm 

2 

16.4Hz 

26.0Hz 

57.1Hz 

68.3Hz 

76.2Hz 

t  t  t 


Ph.  Planquart,  P.  Rambaud 
von  Karman  Institute,  Belgium 


3'^  International  Symposium  on  Integrating  CFD  and  Experiments  in  Aerodynamics 
U.S.  Air  Force  Academy,  USA,  June  20-21,  2007 


Siide7/22 


■  Vortex  shedding  frequencies  have  been  obtained  using  a  CFD  approach. 

■  simulation  of  the  flow  field  for  the 

direction  as  in  the  is  performed. 

■  For  a  ,  the  vortex  shedding  found  numerically  is 

equal  to 


Contours  of  velocity 
magnitude  at  t=1 .1 25  s. 
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u 
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Hz 
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10 
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20 
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67 

0,278 

25 
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Contours  of  veiocity 
magnitude  at  t=1.13  s. 


■  This  has  been  by  the  ,  for  which  vibrations 

related  to  torsion  were  observed  to  begin  at  a  wind  tunnel  velocity  of  10  m/s. 
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■  Ten  different  geometries  have  been  proposed  and  analyzed  (CAE). 


Ph.  Planquart,  P.  Rambaud 
von  Karman  Institute,  Belgium 


3'^  International  Symposium  on  Integrating  CFD  and  Experiments  in  Aerodynamics 
U.S.  Air  Force  Academy,  USA,  June  20-21,  2007 


Siide9/22 


6/7 


Flow  structure  interaction 


■  Natural  of  the  ten  new  blades  were 

■  The  had  to  be  found  (flow  field  around  the  blade). 

■  The  presents  velocity  leading  to  lower  pressure 

fluctuations  and  lower  amplitude  of  vortex  shedding. 


■  The  was  built  and  then  in  the 

■  Very  good  results,  meaning  ,  have  been  observed  at 

low  speed  with  the  proposed  blade  geometry. 
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■  Aerodynamic  of  the  new  for  the 

■  The  presence  of  the  on  Earth,  so  called  the  katabatic  wind 

leading  to  wind  speed  in  excess  of 

■  Transport  of  large  quantities  of  that  can  bury  the  building  very  quickly. 
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Figure  1l  Situation  ciup  «f  thf  facarg  AQ^Hirtic 
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Most  important  aspects  have  to  be  considereo. 


Aerodynamic  design  of  a  new  polar  station 
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Aerodynamic  design  of  a  new  polar  station 


■  The  aerodynamic  design  is  divided  in  two  parts: 

•  A  conceptual  phase 

•  A  detailed  analysis  phase 

■  In  the  ,  the  ability  of  different  buildings  to  cope  with  intense 

drifting  snow  has  been  tested. 

■  This  be  made  with 

■  The  methodology  developed  by  VKI  has  allowed  to  find  experimentally  regions 
of  snow  accumulation  and  snow  erosion. 
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Aerodynamic  design  of  a  new  polar  station 


■  In  the  ,  the  height  and  position  of  the  building  on  the 

ridge  has  been  defined. 

■  on  the  building  has  been  evaluated  for 

of  the  building. 

■  were  also  used  the  pattern  around  the  building. 

The  vortices  behind  the  building  are  responsible  of  the  snow  erosion  and  snow 
accumulation. 

■  have  been  “flow  visualization”. 
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Front  View 


Bottom  View 


Aerodynamic  design  of  a  new  polar  station 


was  built  and  in  the 
in  the  wind  tunnel  have  been 
for  the  of  the  external  panels 


results 


ibution  on  the 
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Aerodynamic  design  of  an  ultra 
streamlined  land  vehicle 


The  is  organized  every  two  years  and  is  a  3000  km  long 

race  for  (www.wsc.org.au/2007/). 

New  for  the  Challenge  the  external  aerodynamic  h  - 

of  the  solar  car.  rf’  ■  P 

Two  phases  for  the  design 


Design  started  by  using  CAD  package, 
the  to  get  aerodynamic  data 


tries  were  then  inserted  in 


Darwin  (0  km) 


Cadmy  Homestead  (2036 km) 


Glendambo  (2442  km) 

>  d  Port  Augusta  (2731 


OT  the  was 

to  define  the  of  the  car 

includes  and 

driver  that  will  be  used  to  power  the  car. _ 


below  the 
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About  5.5  millions  cells 


Aerodynamic  design  of  an  ultra 
streamlined  land  vehicle 


detailed  investigation  of  the  flow  field  around  the  car  (CFD) 
that  were  tested: 

Thickness  and  supporting  the  solar  cells 


is  obtained  in  the 


depending  on 


PRESSmE 


C.OOOeHX 
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Aerodynamic  design  of  an  uitra- 
streamiined  iand  vehicle 


The  following  were 

>  What  is  the 

>  Howto 

>  What  is  the 

>  What  is  the 

>  What  would  be  the 

the  aerodynamic  drag  ? 

>  What  would  be  the 
car  on 


at  the  end  of  the  analysis: 

? 

of  the  to  reduce  the  overall  drag? 

for  the  system  ? 

on  the  drag  and  side  force  ? 
or  dirt  on  the  surface  of  the  solar  cells  on 


to  a  small 


of  the 


had  to  be  performed  because  it  was  the 
to  these  questions. 
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by  CFD  (reduction 


tests  were  performed  on  a 
have  been  performed  in 
(Cd.A)  was 

by  more  than  10%). 
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Aerodynamic  design  of  an  ultra 
streamlined  land  vehicle 


Aerodynamic  design  of  an  ultra- 
streamlined  land  vehicle 

B  WIND  TUNNEL  TESTS 

SIMULATION 

1  Wind  tunnel  tests  First  estimation  of  the 

1  with  two  different  models  aerodynamic  coefficients 

1  Finalization  of  the  design 

Parametric  tests  to  create 
a  database  of  aerodynamic 
coefficients 
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Conclusions 


■  can  lead  to  a 

■  wind  tunnel  and  CFD  specialists 

throughout  the  studies. 

■  The  have  also  shown  that 
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